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Abstract 

Motivated by a recent discussion about the role of flat directions, a typical feature of 
supersymmetric models, in the process of particle production in the early universe a consis- 
tent model of inflation and preheating in supergravity with MSSM fields has been built. It 
is based on a model proposed by M. Kawasaki, M. Yamaguchi and T. Yanagida. In the in- 
flationary stage, the flat directions acquire large vacuum expectation values (VEVs) without 
spoiling the background of slow-roll, high-scale inflation consistent with the latest WMAP5 
observational data. In the stage of particle production, naturally following inflation, the role 
of flat direction large VEVs depends strongly on effects connected with the supergravity 
framework and non-renormalizable terms in the superpotential, which have been neglected 
so far in the literature. Such effects turn out to be very important, changing the previous 
picture of preheating in the presence of large flat direction VEVs by allowing for efficient 
preheating from the inflaton. 



1 Introduction 

Inflation was introduced as a natural and simple way of solving the problems of classical cosmol- 
ogy - the initial conditions problem (or the flatness and horizon problems) and the explanation 

of the origin of primordial density fluctuations [UEJEJIHEJEIEJE]. ^ 6 6aS * 6S ^ wa ^ ^° °^^ am 
inflation is by introducing a single scalar inflaton field with a slowly evolving vacuum expectation 
value [2j [3]. In order to obtain a proper period of Big Bang Nucleosynthesis however, one has to 
end inflation by particle production. The process of reheating must connect the inflaton sector 
to the observable sector [91 fTTJl [TT] . 

In order to properly describe inflation and particle production one has to consider the underlying 
theory of particles and interactions. Supersymmetry is one of the most promising extensions of 
the Standard Model (SM) [I2j [TH1 HH HH HH [TT] , and it has triggered a search for supersymmet- 
ric models of inflation and reheating. One of the typical features of supersymmetric extensions 
of the SM is the presence of flat directions [TH] - directions in field space, along which the scalar 
potential identically vanishes in the limit of unbroken global supersymmetry. Due to large quan- 
tum fluctuations or the classical evolution of fields during inflation flat directions can easily 
acquire large VEVs [T5J[20]. Therefore, there is a natural question about the role of such large 
VEVs in the process of particle production. 

It was postulated in ref. [2Tj, that large flat direction VEVs influence the process of particle pro- 
duction by blocking preheating from the inflaton - the phase of rapid, non-perturbative inflaton 
decay. In ref. [21J a simple toy model was proposed 

V D Aip 2 X 2 + Brrnpx 2 + Ca\ 2 , (1) 
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where tp is the inflaton field, a parameterizes the flat direction and x represents the inflaton 
decay products (in this model a direction in Higgs fields has been considered). Then, after mode 
decomposition of the field Xi the energy of the mode with momentum k is given by: 

co 2 k = k 2 + 2A(<p} 2 + 2Bm(ip) + 2C (a) 2 . (2) 

In general, non-adiabatic production of particles Xk is efficient only when Wfc changes non- 
adiabatically 



UJ 2 



> 1 <-» preheating, (3) 



where the adiabaticity parameter r is introduced. During classical preheating ujk is dominated 
by the inflaton VEV and changes non-adiabatically due to inflaton oscillations. In the presence 
of flat directions however, u>k could be dominated by the large VEV of the flat direction. If this 
VEV changes very slowly in comparison with the evolution of the inflaton VEV, non-perturbative 
production of x particles is effectively blocked. 

However, as was pointed out in ref. [22], blocking of preheating from the inflaton does not occur 
when non-perturbative production of particles from the flat direction itself is possible. Then 
the initially large VEV of the flat direction decreases rapidly, unblocking preheating from the 
inflaton. In ref. [22] a method was introduced of calculating the amount of particles produced 
non-perturbatively from the flat direction, due to non-adiabatic changes of the mass matrix 
eigenvectors and eigenvalues related to quantum fluctuations around the flat direction. This led 
to a discussion (see refs [231 EH ESI EE]) about whether non-perturbative decay of flat directions 
and preheating from the inflaton is possible. The discussion was based on some general properties 
of flat directions in a global supersymmetry framework. It did not consider any specific model 
of inflation and did not propose any model of acquiring large VEVs by flat directions. Therefore 
it was difficult to study the whole issue and determine how the large VEVs of flat directions 
develop and evolve, and how they impact the process of inflation and particle production. 
The goal of our work is to construct a consistent model of inflation and particle production 
in a realistic supersymmetric extension of the Standard Model, and consider in this specific 
model the behavior of MSSM flat directions. Therefore a realistic chaotic inflation model with 
two representative flat directions is constructed. In order to be able to predict the evolution 
of flat direction VEVs it was decided to study the production of large flat direction VEVs by 
classical evolution during inflation, and therefore a potential for the flat direction is required. 
Following refs [19, 20J we adopt the supergravity framework with a non-minimal Kahler potential, 
which results in a potential for the flat direction with a time-evolving minimum at large VEVs 
during inflation. It also enables us to calculate the previously neglected influence of supergravity 
corrections for flat direction evolution during inflaton oscillations. We also consider the impact of 
existence of non-renormalizable terms, which has also been neglected so far. We find that these 
effects strongly influence the process of particle production by introducing efficient channels of 
non-perturbative particle production both from the flat direction and the inflaton. As a result 
the originally large flat direction VEVs are diminished, preheating from the inflaton is allowed 
and the energy density of the Universe is dominated by the inflaton decay products. 



2 Building the model 

The model considered in this paper, after neglecting all fields except for the inflaton, reduces 
to the simplest chaotic inflation model with the inflaton potential V = m 2 (p 2 /2. This property 
together with inflaton domination provides appropriate slow-roll inflation and a value of spectral 
index which is in agreement with the WMAP5 data [27]. Obtaining such a property in a model 
which is based on supergravity is not straightforward due to complicated supergravitational F- 
terms. A solution to this problem (the so called 77-problem) was proposed by [28] and is used in 
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this paper. According to the solution, except for the chiral inflaton superfield $ and the MSSM 
superfields, the model contains one additional chiral superfield X. 

Further consideration is restricted to scalar fields and the same symbol is used to denote both 
the chiral superfield and its complex scalar component. The following decomposition in real 
fields is used 

$ = (7/ + i<p)/y/2, (4) 
X = xe ip . (5) 

The field <p plays the role of the inflaton. 

We follow [28] in constructing the Kahler Potential K and take 

Kd -($ + $*) 2 + XX*. (6) 

The formula (JSJ) for the Kahler potential was obtained by [28] as follows. The first step was to 
introduce a Nambu-Goldstone-like shift symmetry of the inflaton a 

$ -► $ + iCM, (7) 

where C is a dimensionless real parameter. A Kahler potential which is invariant under this 
symmetry and the additional U{1)r x Z2 symmetry must have the general form [28] 

= K[($ + $*) 2 ,XX*}. (8) 

The formula J6]) is just the lowest order term in the general expansion of the formula It 
can be easily seen that a theory with an exact Nambu-Goldstone-like shift symmetry has no 
potential for the inflaton (p. Therefore this symmetry has to be broken, but not in the Kahler 
potential in order to avoid the ry-problem. Following [28] . we introduce a shift symmetry breaking 
term in the superpotential W 

W D ml§. (9) 

It gives mass m to the inflaton ip. 

There are two MSSM-flat directions considered in this paper: a direction H u Hd in Higgs fields 
(only D-flat) and a direction Uidjd^ in squark fields, where indexes i, j and k are some family 
indexes (k / j). Let x t> e the complex scalar field that parametrize the H u Hd direction 

Let a be the complex scalar field which parametrizes the Uidjdk direction 



- XI - 
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d] = d{ = —a. (II) 



V3 



In equation (fTTj) (3 ^ 7 / 5 / (3 are fixed color indexes. The components of fields «j, dj and d^ 
with other color indexes are equal to zero. It is convenient to decompose the complex fields \ 
and a into real fields in the following way 

X = ce lK , (12) 
a = pe ia . (13) 



The full Kahler potential is 



K = i($ + <S>*f + XX* + Kmssm + K NM , (14) 
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where Kmssm is a standard minimal Kahler potential and Knm is a non-minimal part of the 
form 

K NM = jfcXX*- 

■(H+^ + H+Hd + u+Ui + d+dj+d+dk). V ' 

Here Mpi is the Planck mass and a is a dimensionless parameter. The existence of couplings like 
Knm is guaranteed in the presence of Yukawa couplings, since they are necessary counterterms 
for operators generated by loop diagrams [1J2 [29j [30]. Terms in Knm cause the existence of 
minima in the scalar potential for both flat directions, which are of the order of Mp\. Therefore 
flat directions can naturally aquire large VEVs of the order of Mp\ by falling into these minima. 
For the superpotential we take 

W = ml$ + 2hXH u ■ H d + Wmssm + W NR . (16) 

The term 2hXH u ■ Hd is one of only two possible renormalizable couplings between the field 
X and MSUGRA fields. The second one is 2h'XH u • L. These two couplings cannot coexist in 
the model unless R-parity is broken. We assume that the inflaton does not couple to MSSM 
fields in the Kahler potential and in the superpotential to avoid strong deviations from the 
slow-roll inflation regime with inflaton domination. The term Wmssm is the standard MSSM 
superpotential, given by 

Wmssm = »H U ■ H d + X l ™H u ■ Q x u m + X l ^H d ■ Qid m . (17) 

The last term Wnr is a non-renormalizable part of the superpotential and it has the following 
form 

W NR = ^(H U - H d ) 2 + (u t djd k v R ) . (18) 

Mpi Mpi 

Here A x and X a are dimensionless constants and up is a right-handed neutrino of any given 
generation. Two non-renormalizable terms contained in Wnr are the only terms of 4th order in 
the fields, which may be relevant for the evolution of VEVs of the two chosen flat directions. 
The term 3\^3X' a (uiUjdkVp) /Mp\ gives no contribution to the scalar potential, unless one ad- 
ditionally considers VEVs of some other flat directions. 

Terms in Wnr modify minima for flat directions, which are shifted away from zero in the scalar 
potential due to terms in Knm- The minima are no longer close to Mp\ when coupling constants 
X x and X a are sufficiently large. Moreover, they evolve in time during inflation until they reach 
their final position at zero at the end of the inflaton oscillations. 

The soft SUSY-breaking terms in the scalar potential have negligible effects on our results. 
Our initial conditions are set for the time which corresponds to about 100 e-folds before the 
end of inflation, since only this period is essential for preheating. Initial values for the real fields 
ip, rj, x, p and c do not require fine tuning. For the inflaton field ip the only condition which 
needs to be satisfied is cpo > Mp\. It ensures that the number of e-folds is greater than 75, which 
is needed for inflation to solve the horizon problem and the flatness problem |8|. We took the 
initial value (po = 4Mpi to consider only the last 100 e-folds before the end of inflation. Fields 
rj, x, p and c should be initially smaller than Mpi because they are present in Kahler potential, 

especially in the exponential factor exp(yp-) in the F-terms. Therefore initially large values of 

m pi 

any of those fields should decrease rapidly and the VEVs of the fields t], x, p and c should stay 
confined naturally below the Planck scale during inflation. In particular the VEV of field r\ falls 
to zero very quickly and, as we have checked numerically, does not have any noticable influence 
on the evolution of other fields. Therefore in our final calculations we simply put rjo = 0. The 
evolution of the field x has also almost no influence on the evolution of other fields, so we set 
for it a quite arbitrary initial value xq = O.OIMp/. Fields p and c which are absolute values for 
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complex flat direction fields a and x are initially taken to be of the order of the Hubble param- 
eter H, which is also the order of initial quantum fluctuations for those fields during inflation. 

In the slow-roll regime with inflaton domination, we have Hq 87rm 2 (p\ / 6Mj, ; . After setting 

m = W~ 6 Mpi as in the simplest model of slow-roll inflation consistent with WMAP data we 
get H Q « 8 • 1(T 6 Mpi. Therefore we take p = c = 10~ 5 M P/ . 

3 Classical evolution of fields 

In order to study the process of acquiring large VEVs by flat directions one has to consider the 
classical evolution of fields during inflation. A classical description is possible due to the slow-roll 
character of the evolution. At the end of inflation, excitations around VEVs of the fields will 
be considered in order to determine the impact of large flat direction VEVs on the process of 
particle production. 

The classical evolution is determined by the equations of motion derived from the supergravity 
Lagrangian once the initial conditions have been set. The inflaton equation of motion has the 
simple form 

(p + ZHip + V, v = 0. (19) 

During inflaton domination the main contribution to the scalar potential V is of the form 
1/2 m?Lp 2 , which provides a standard chaotic inflation background. Due to the shift symme- 
try the Kahler potential does not depend on ip, and so the supergravity coefficient e K does not 
contain the inflaton field. This solves the r/-problem and allows for high-scale inflation. Since 
only the last 80-100 e-folds of inflation have any observable consequences for the evolution of 
the Universe, the initial value of the inflaton VEV was chosen in a way that allows the study 
of this period of inflation. Due to the absence of any coupling of the inflaton with the H u Hd 
and udd directions in the Kahler potential (because of the shift symmetry) the evolution of the 
inflaton is largely independent of the evolution of the flat directions. Because of this property 
the process of acquiring large VEVs by udd or H u Hd directions will not spoil inflation. 
The equations of motion for X, udd and Higgs fields are more complicated due to the non- 
minimal form of the Kahler potential for these fields. One can expect however that the field 
VEVs will evolve toward a minimum of the scalar potential. The scalar potential in supergravity 
is complicated as well. However one can observe that all the field VEVs except for the inflaton 
VEV are naturaly confined below the Planck scale due to the e K factor in the scalar potential. 
This factor becomes dominant in the scalar potential at the Planck scale. Therefore the scalar 
potential for all fields except the inflaton rises steeply at the Planck scale as exp(field 2 ). As a 
result one can expect x, p, c to be less than unity and expand the scalar potential in these fields. 
For field X = xe % ^ the scalar potential has a minimum at zero. The term quadratic in field x 
is given by the following approximate expression (assuming inflaton domination and neglecting 
complex fields phases) 

V D x 2 (m 2 + (f(a)f? + / (a)c 2 + O [p\\ p\ c 3 ])) , (20) 

where for simplicity we have set Mpi = 1 and /(a) is positive for a > 0. The evolution of x is 
naturally confined to low VEVs due to the supergravity term e K which exponentially steepens 
the potential for field x at the Planck scale. The evolution of all other fields of the model is 
independent of the choice of initial conditions for x and on the evolution of this field. The only 
role of field X in the model is providing appropriate scalar potential for the inflaton and both 
flat directions. 

In the global MSSM without corrections coming from supergravity or non-renormalizable terms 
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the scalar potential is independent of the flat direction a = pe la . Adding the corrections men- 
tioned above creates a potential for the flat direction. A term quadratic in a in the scalar potential 
is created by supergravity effects and is sensitive to any non-minimal couplings in the Kahler 
potential. In the model presented in this paper the quadratic term mentioned above takes during 
inflaton domination (neglecting the complex fields phases) the following approximate form 

V D p 2 ^f- (1 - a + f{a)c 2 + f(a)x 2 + O [x 2 c 2 , x 3 , c 3 ] ) . (21) 

In the equation above the parameter a describes the influence of the non-minimal coupling 
in the Kahler potential. In supergravity with a minimal Kahler potential the scalar potential 
has a global minimum at zero for the flat direction. In this case acquiring large flat direction 
VEVs due to classical evolution is impossible. The non-minimal coupling enables us to shift 
the minimum toward larger VEVs by an appropriate choice of a > 1. Then the coefficient of 
the term quadratic in p in the scalar potential becomes negative for x, c < 1. For the purpose 
of numerical calculations we set a = 5. The exact location of the minimum is determined by 
higher-order terms in p, which come from both supergravity and non-renormalizable terms. 
Supergravity alone stabilizes the minimum around Mp\ due to the coefficient e . The presence 
of a non-renormalizable term scaled by X a shifts this minimum toward lower VEVs, changing the 
predicted flat direction VEVs at the end of inflation toward lower values. In order to study the 
effect of large flat direction VEVs on the process of particle production the value of X a should 
not exceed unity. 

Supergravity corrections and non-renormalizable terms have a similar effect on the potential for 
the H u Ha direction. The interplay between A x and X a determines the differences between the 
evolution of the two flat directions (H u Hd and udd) under consideration. We study two specific 
scenarios 

1. If A a <C A x then the udd direction VEV becomes large and the HuH^ direction VEV 
drops to zero during inflation, which corresponds to the scenario described in ref. [21] in 
the limit of global supersymmetry without non-renormalizable terms. This case will enable 
us to study the predictions of ref. [21] in a specific scenario and determine the impact of 
supergravity corrections and non-renormalizable terms, which has not been considered so 
far. 

2. If A a ~ A x then both directions can acquire large VEVs during inflation and the impact of 
non-zero VEVs of the inflaton classical decay products on the process of particle production 
and blocking of preheating by flat direction large VEVs can be studied. 

By "large VEV" we mean a vacuum expectation value of the order 10~ 4 — IMpi, which is large in 
comparison with the Hubble parameter H ~ 1CT 7 Mpi at the beginning of inflaton oscillations. 
Since the acquired value of the uud or H u Hd direction VEV is determined by the value of A a 
and A x respectfully, in order to create large VEVs of those directions during inflation X a and A x 
should not exceed unity. For the purpose of numerical calculations two sets of A-parameters are 
considered 

1. A a = 1CT 7 and A x = 1 

2. X a = 1 and A x = 1. 

In order to check if classical evolution can lead to large flat direction VEVs during inflation 
we have chosen small initial VEVs for both directions, which correspond to the average size of 
quantum fluctuations typical for the considered period of inflation (5a, 5c ~ H). In order to 
obtain numerical predictions in specific scenarios one has still to fix two free parameters. The 
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choice of the inflaton mass m ~ 10~ 6 is natural because it implies a spectral index of energy 
density fluctuations consistent with WMAP observations. Following the arguments of ref. [31] 
we choose the inflaton coupling parameter h to be small and of the order h ~ 10 -5 . 
The classical evolution of fields obtained numerically in the two scenarios mentioned above is 
presented below. It turns out that in both cases the evolution of field rj is irrelevant - the VEV 
of the 77 field decreases rapidly to zero and does not influence the further evolution of the other 
fields. Henceforth we simplify our calculations by setting r\ = 0. 

3.1 X a = lCT 7 and A x = 1 

Preliminary numerical calculations show that in this case some of remaining 7 real VEVs (after 
setting 77 = 0) can be also neglected. First of all the scalar potential does not depend on the 
phase a and any initial velocity of this phase will fall to zero due to the Hubble friction term. 
Therefore the flat direction udd can be effectively described by only its absolute value p. The 
VEV of the field c, the absolute value of Higgs flat direction, after an initial increase, falls to 
zero quite rapidly. Thereafter the scalar potential ceases to depend on f3 and k. 
Henceforth we present more accurate calculations only for 4 real fields: </?, x, c and p. All three 
phases have been set to zero. In this scenario the final 100 e-folds of inflation are studied. 
The inflaton VEV evolves according to standard chaotic inflation. It smoothly decreases during 
inflation and starts to oscillate after its end, as can be seen in figures JT]) and (J23) - The field x 



Figure 1: Evolution of the inflaton field <p dur- 
ing inflation. Values on vertical axes are ex- 
pressed in Planck masses and time on hori- 
zontal axes is expressed in the approximate 
number of e-folds of inflation. 




Figure 2: Inflaton oscillations at the end of in- 
flation. Values on vertical axes are expressed 
in Planck masses and time on horizontal axes 
is expressed in Planck times. 



behaves almost identically to the inflaton field (p. The main difference is that the VEV of the field 
x is much smaller than the VEV of the inflaton. This is illustrated in figures j3|) and (|4}. The 
VEV of the field c initially rises due to the influence of the non-minimal Kahler coupling which 
shifts the minimum of the scalar potential for this direction away from zero. After the initial 
rise this VEV starts decreasing and drops to zero rapidly during the first 35 e-folds of inflation. 
This effect is shown in figure j5|) and is caused by the relatively large value of A x (A x = 1) with 
respect to A a , which favors the creation of large udd direction VEV. Due to the Yukawa coupling 
between udd and H u H<i directions a large udd direction VEV induces an effective mass for the 
H u Hd direction. This effect shifts the minimum of the scalar potential for the H u Hd direction 
to zero leading to the rapid decrease of the H u Iid VEV during inflation. 

The VEV of the field p rises during inflation to a value which is close to Planck mass due to the 
small value of \ a (Fig. J6])). At the end of inflation it starts to slowly decrease, as illustrated in 
Fig. (J7|. The evolution of the spectral index in the crucial period 60-50 e-folds before the end 
of inflation can be calculated in the slow-roll regime and is shown in figure JH)). The value of the 
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Figure 3: Evolution of the field x during in- 
flation. Values on vertical axes are expressed 
in Planck masses and time on horizontal axes 
is expressed in the approximate number of e- 
folds of inflation. 



Figure 4: Oscillations of the field x at the 
end of inflation. Values on vertical axes are 
expressed in Planck masses and time on hor- 
izontal axes is expressed in Planck times. 
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Figure 5: Evolution of the field c during during inflation. Values on vertical axes are expressed 
in Planck masses and time on horizontal axes is expressed in the approximate number of e-folds 
of inflation. 




Figure 6: Evolution of the field p during in- 
flation. Values on vertical axes are expressed 
in Planck masses and time on horizontal axes 
is expressed in the approximate number of e- 
folds of inflation. 
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Figure 7: Evolution of the field p at the end 
of inflation. Values on vertical axes are ex- 
pressed in Planck masses and time on hori- 
zontal axes is expressed in Planck times. 



spectral index 50 e-folds before the end of inflation is in agreement with the value derived from 
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Figure 8: Evolution of the spectral index n$ between 100 and 40 e-folds before the end of 
inflation. Black lines marks the time of 50 e-folds before the end of inflation and corresponding 
value of the spectral index. Time on horizontal axes is expressed in Planck times. 

the WMAP5 data n s = 0.960lg;gJ| [21]. 
3.2 Aq, — 1 and A x = 1 

The inflaton VEV evolves according to standard chaotic inflation as in the previous case. Its 
evolution includes the slow-roll period which naturally ends with inflaton oscillations. Due to 
the dependence of the scalar potential on the phases of fields X and x> both of these fields evolve 
non-trivially in the complex plane. The evolution of the absolute value of field X is similar to 
the previous scenario and mimics the behavior of the inflaton. Figure j9]) shows the evolution 
of field X in the complex plane. Since X a ~ A x now the absolute values of fields corresponding 



Figure 9: Evolution of field X on the complex plane during inflaton oscillations. 

to both flat directions evolve similarly. Numerical calculations show that the absolute values of 
both fields grow during inflation, achieving a maximum of the order of 10 _3 Mp; and then start 
to decrease at the end of inflation due to the time-evolution of the minimum of their potentials. 
Because of the dynamics of the phase of field x the effective mass of the field is slightly different 
to that of the field a. As a result field x begins oscillating earlier. Figures (fTOjl and ([Pil show the 
evolution of the absolute value of x, while figure (TT21) shows the evolution of x m the complex 
plane. The evolution of the udd flat direction is effectively one-dimensional because the scalar 
potential does not depend on the phase of this field and the Hubble friction quickly suppresses 
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Figure 10: Evolution of the absolute value of 
the Higgs field direction during inflation. Val- 
ues on vertical axes are expressed in Planck 
masses and time on horizontal axes is ex- 
pressed in the approximate number of e-folds 
of inflation. 
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Figure 11: Evolution of the absolute value of 
the Higgs field direction during inflaton oscil- 
lations. Values on vertical axes are expressed 
in Planck masses and time on horizontal axes 
is expressed in Planck times. 
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Figure 12: Evolution of the Higgs field direction in the complex plane during inflaton oscillations 



any initial phase dynamics. Figures (TT3l) and JHJ show the evolution of the absolute value of a. 
Using the numerically found evolution the Hubble parameter and slow-roll parameters can be 
calculated. They fulfill the slow-roll conditions. One can then obtain the spectral index in the 
slow-roll approximation. The spectral index evaluated at the time of about 50-60 e-folds before 
the end of inflation is consistent with the value of the spectral index derived from the WMAP5 
observational data |27j . 

4 Excitations around VEVs 

In this chapter we introduce excitations around all MSUGRA fields, which are related to the 
flat directions under consideration. By "related" we mean that either they have large VEVs, 
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Figure 13: Evolution of the absolute value of 
the udd flat direction during inflation. Val- 
ues on vertical axes are expressed in Planck 
masses and time on horizontal axes is ex- 
pressed in the approximate number of e-folds 
of inflation. 
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Evolution of the absolute value of 



the udd flat direction during inflaton oscilla- 
tions. Values on vertical axes are expressed in 
Planck masses and time on horizontal axes is 
expressed in Planck times. 



parametrized by flat directions, or they are parts of multiplets in which other fields have such 
VEVs. In the first case we parametrize those fields in the following way 

VEV\ h \ d„,(vEV >+7 ^) 



Fl = {^r + f2) e (22) 

In the second case the parametrization has the form 

F2 = ~fa + i$ 4 )e i - ar <><y EV ). (23) 

In formulas (|2"2"|) and (|2"3"1) £i, £2, £3 and £4 are real excitations, whereas VEV denotes x or 
a. Moreover n = 2 for Higgs doublets and n = 3 for squark triplets. If the VEV of a flat 
direction drops to zero, the parametrization of excitations around the related MSUGRA fields 
is straightforward 

F 3 = -J= (6x + id 2 ) . (24) 
Excitations around fields $ and X are not considered. 

We initially have 26-dimensional space of real excitations. Some are Goldstone bosons related 
to the gauge group generators which are broken by the VEVs of flat direction(s). Goldstone 
bosons can be eliminated via the Higgs mechanism. Then, expanding the Lagrangian density in 
the remaining excitations 

L = L + L x + L 2 + . . . (25) 

Here L n is the part of Lagrangian density which is of n-th order in the excitations. In particular 
Lq is the classical and homogeneous limit of the Lagrangian density, used previously to obtain 
the evolution of VEVs. Terms L n for n > 2 are neglected from now on. L\ can be set to zero after 
using partial integration of the action and the classical equations of motion. The remaining part 
L 2 is used to describe particle production on the classical, homogeneous background including 
the dynamics of VEVs and the scale factor. The kinetic part of L 2 is quite complicated due 
to the non-minimal part in the Kahler Potential K^m- However, as was mentioned in section 
3, excitations are considered only at the end of inflation and later. In this period Knm can 
be neglected since the multiplicative factor included in equation (TT5j) . is very small. 
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This approximation is adopted only for kinetic terms, keeping the full SUGRA scalar potential 
(we checked numerically that using minimal kinetic terms during inflaton oscillations does not 
change the classical evolution of the fields, while changing the potential does alter the evolution). 
Then, L 2 has the following general form 

L 2 = \d»E T d»E - \e t M 2 E - E T UE. (26) 

Here H is the vector which contains all excitations and M 2 is the mass matrix. There is also a 
matrix U, which mixes excitations with their derivatives. Elements of the matrices M 2 and U 
can in general be functions of all 8 real VEVs and the matrix U can also contain derivatives of 
those VEVs. There are two steps needed to transform L 2 into more convenient form: 

1. Integrating the term E T UE by parts in the action in order to replace the matrix U with 
an antisymmetric matrix U. This procedure gives new contributions to the mass matrix. 
Its new form will be denoted by M 2 

L 2 = \d»E T d»E - \E T M 2 E - E T UE. (27) 

2. Since U is antisymmetric, one can find an orthogonal matrix A such that 

U = A T A. (28) 
Defining H' = AE, the U matrix can be eliminated 

L 2 = 1^H*>3' - l -E' T M' 2 E\ (29) 

where M' 2 = A (M 2 - U 2 ^j A T . 

Similar transformations of L 2 were presented in [25]. The matrix M' 2 can be diagonalized 

M' 2 = CMjC T , (30) 

where the matrix Mj is diagonal and the matrix C is orthogonal. The situation is more general 
than in [25], because both these matrices are functions of VEVs, so both of them depend on 
time. We have time dependent eigenvectors, but also time dependent eigenvalues of matrix M' 2 
and both these time dependences influence particle production. 

To calculate particle production we quantize excitations in curved space-time according to [3^] . 
Quantum excitations are minimally coupled to gravity (they don't have any couplings to the 
Ricci scalar) and have effective squared masses which may differ significantly from the eigenval- 
ues of M' 2 , according to the following formula 

ml m = m| - 2H 2 - H. (31) 

Here m| is the ordinary squared mass of an excitation £ (an eigenvalue of M' 2 ), m 2 ^^ is the 
effective squared mass of this excitation and H is the Hubble parameter. 

Two sets of excitation modes are used so that the corresponding vacua minimize Hamiltonian 
in two particular moments of time. The set of in modes minimize the Hamiltonian at to, the 
beginning of the considered particle creation period, and the set of out modes minimize the 
Hamiltonian at tj, the end of this period. After finding numerically the evolution of in modes 
between to and t\ the Bogolyubov coefficients method [32j [33j [34] is used to obtain the energy 
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density of produced particles. It is worth noting that the above method allows us to describe 
particle production in the fully non-perturbative regime. We do not need space-time to be 
initially close to de Sitter and we do not use adiabatic modes. Non-adiabatic particle production 
is efficient when the adiabaticity condition j3]) is broken, [21]. Therefore in each of the cases, 
we choose a period of time for our calculations that corresponds to the broken adiabaticity 
condition. 



4.1 X a = 10" 7 and A 



x 



In this case the only fiat direction which still has a VEV at the end of inflation is the udd 
direction. Therefore excitations which correspond to Higgs doublets are described by formula 
(|24l) . Other excitations, which correspond to squark triplets, are described by formulae (|22j) and 
(|23l) . The VEV of the field p breaks the gauge symmetry 



SU(3)c x U(1) Y - U(1) P . (32) 
Here U(l)p is parametrized by a single generator P, which is defined in the following way 

P = Y-\j 3 -^J s , (33) 

where Y is the weak hypercharge, while J3 and Jg are two SU (3)c generators, which can be 
represented as two diagonal Gel-Mann matrices. The SU(2)l symmetry remains unbroken. There 
are 8 broken generators, related to 8 Goldstone bosons among 26 initial excitations. Hence there 
are 18 physical degrees of freedom in the excitation space. After eliminating Goldstone bosons, 
we use the unitary gauge to parametrize this 18-dimensional space. The excitations around fields 
in the Higgs doublets are 

H ul = ■j^{5 l + i5 2 ), 
H u 2 = ^(63 + 164), 

H d2 = -^(Sr + iSs). 

Excitations around fields in squark triplets in the unitary gauge (after eliminating Goldstone 
bosons) take the form0 

Ui 2 = \ (& + *&) e la , 
u i3 = 5(60 + e i0 ", 
dj 1 = \ (& - «&) e i<J , 
^2=f^ + % N )e< CT+ 7fe), (35) 



4 1 = 1 (60 - ibx) e i<T , 

42 = 1 (63 - %u) e ia , 

*3 = (^ + H)e<' + *). 

We followed the procedure described earlier in this section and found that the U matrix van- 
ishes in this case, because there is no phase dynamics, so M' 2 = M 2 . The mass matrix of the 
excitations related to Higgs fields and udd fields has a block diagonal form 

M' 2 - ( M 8x8 [ H uH d ] \ . . 

V M 2 0xl0 [udd]) ■ (36) 



1 In section 3.1. a has been set to zero 
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The 8 x 8-dimensional block of Higgs-related excitations has 8 eigenvalues, which (under the 
simplifying assumption that the Yukawa matrix is diagonal in flavor and X u = Xd = Y) are 
degenerate. Four of them have the following approximate form 



~- (2V2h + (a - 1) rrup^j 

+ + V2 (a - 1) hrrnp + (l - 2a + 2a 2 ) p 2 + Q ^ ^,2^ p3 j ^ (3?) 

while the other four are given by the following approximate expression: 

m 2 ~ (-2\/2/i + (a - 1) m<p) 

/y2 2 2\ 

+ f — - V2 (a - 1) hrmp + (l - 2a + 2a 2 ) — ^- J p 2 + O [x 2 , x 2 p 2 , p 3 ] , (38) 

where we have neglected the soft masses and the Higgs mass parameter p. In the equations above 
it can be seen that there are, as expected, two main contributions to the mass eigenvalues. One 
comes from the standard 3- and 4-linear interaction terms with the inflaton. The other is related 
to the flat direction and the dominant contribution comes from the Yukawa interaction between 
Higgs fields and the udd flat direction fields. The first term in each parenthesis, scaled by h or Y, 
is related to the assumed couplings in the superpotential and would survive even in the absence 
of gravitational effects. All the additional terms come from supergravity and the parameter a 
scales the influence of the non-minimal coupling in the Kahler potential. Due to the Yukawa 
coupling the influence of the udd flat direction VEV on the mass eigenvalues related to Higgs 
fields dominates. Therefore these eigenvalues are large and evolve slowly in time, so preheating 
from the inflaton into particles related to those eigenvalues is initially blocked. 
The 10 x 10-dimensional block of the mass matrix, which is related to udd excitations, is also 
block-diagonal. It has one 3-dimensional block and seven 1-dimensional blocks. 

/M 2 xl [phase] 0\ 

Ml x3 [flat] 

A*i xl [l] 

M 2 xl [2] 

A*i xl [3] 

M 2 xl [A] 

M 2 xl [5] 

\0 00000 M 2 xl [6}/ 

(39) 

Eigenvalues M 2 [1] - M 2 [6] correspond to combinations of excitations around fields with VEV = 
belonging to udd flat direction multiplets. They are all heavy because they are related to 
Higgs particles coming from broken non-diagonal generators of SU(3). For example (under the 
simplifying assumption that all the gauge couplings are equal, gi = g) 

M 2 [1] « -Hff (a - 1) + + (1 - 2a + 2a 2 ) + O [x 2 ]) p 2 + O [x 2 , p 3 ] . (40) 

It can be observed that the contribution to this eigenvalue coming from the inflaton VEV is a 
supergravity effect. This is true for all the eigenvalues of M 2 [udd] - in the global SUSY case 
udd fields are not coupled to the inflaton. The dominant contribution to M 2 [1] comes from 
D-terms, is scaled by the gauge coupling g and is proportional to the VEV 2 of the flat direction. 
Therefore all the eigenvalues M 2 [1] - M 2 [6] related to broken non-diagonal generators are heavy 



M [udd] = 
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and evolve slowly in time. At the beginning of inflation there is no non-perturbative particle 
production of particles corresponding to those eigenvalues. The block M 2 [flat] corresponds 
to excitations £7, £12 and £15 around the absolute value of the VEV of flat direction udd. It 
can be diagonalized easily, giving two heavy eigenvalues and one light eigenvalue. Both heavy 
eigenvalues are dominated by terms ~ g 2 p 2 and are Higgs particles corresponding to two broken 
diagonal generators of SU(3). As in previous cases preheating into those particles is blocked. 
The light eigenvalue m 2 corresponds to the combination (£7 + £12 + £15) /Vs of excitations and 
is given approximately by 

22 22 

m | ~ J<J^_ ( fl _ ^ + 3 ^ _ 2a + 2a 2) p 2 + Q ^,2^ ^3] _ ^ 

Clearly this eigenvalue is dominated by supergravity effects (all terms written explicitly in eq. 
(|4~T1) are induced by supergravity) - the only contribution which would remain in global SUSY 
comes from the non-renormalizable term and is of the order ~ p 4 . In global SUSY without non- 
renormalizable terms this eigenvalue would be equal to zero. This is easy to understand since 
in global SUSY without non-renormalizable terms the scalar potential does not depend on the 
flat direction. Then there exist two massless eigenvalues corresponding to excitations around the 
complex field a parameterizing the flat direction. When the scalar potential depends on the flat 
direction (due to non-renormalizable terms or supergravity effects), those two states gain mass. 
The eigenvalue m| corresponds to one of such states. The mass eigenvalues corresponding to 
excitations around a remain naturally light, because any contributions to them originate from 
non-renormalizable or supergravity terms and are diminished by the Planck scale. Moreover, as 
can be seen from (JHJ, the influence of the inflaton VEV on these eigenvalues can easily become 
dominant. Numerical calculations show that the time evolution of this eigenvalue allows for both 
efficient preheating from the inflaton and efficient non-perturbative particle production from the 
flat direction into excitations around the flat direction. Without supergravity corrections or non- 
renormalizable terms this effect would not be possible - the discussed mass matrix eigenvalue 
would be equal to zero. Adding a non-renormalizable term for the flat direction in global SUSY 
would make the eigenvalue non-zero and equal 15X 2 p 4 . Such an eigenvalue is light and depen- 
dent only on the flat direction VEV. This would lead to non-perturbative particle production 
from the flat direction due to the time evolution of this eigenvalue. Adding supergravity effects 
couples the flat direction to the inflaton, allowing for non-perturbative particle production from 
the inflaton as well. 

The second light eigenvalue of the mass matrix, which would be equal to zero without super- 
gravity corrections or non-renormalizable terms, is M 2 [phase] and is related to the excitation £2 
around the a phas^E The squared mass M 2 [phase] = m| 2 (an element and eigenvalue of matrix 
M' 2 ) has the following form 



2 p + 3Hp V, p 
m * = — = -2p~ (42) 



« (1 - a) ^f- + (1 - 2a + 2a 2 ) ^f- p 2 + O [x 2 , x 2 p 2 , p 3 ] , (43) 

where H is the Hubble parameter. Due to the simple form of this mass eigenvalue, which comes 
only from the contribution of the kinetic terms to the mass matrix, one can write this eigenvalue 
explicitly. The mass eigenvalue m| 2 is the smallest of the eigenvalues of the matrix M' 2 . It can be 
observed that the eigenvalue is dominated by supergravity terms. The non-renormalizable term 
gives a contribution 2>X 2 p p A . In global SUSY without non-renormalizable terms this eigenvalue 
would be identically equal to zero. Evolution of the adiabaticity parameter J3|) corresponding 
to m^ 2 at the end of inflation is shown in Fig. (fT5l) . It can be seen that the behaviour of the 



initially, there are 3 excitations around <r, but two of them are Goldstone bosons. 
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adiabaticity parameter becomes quickly dominated by the influence of the inflaton oscillating 
VEV. This leads to the breaking of the adiabaticity condition. 

To calculate production of particles one has to derive an effective mass squared according to 
equation (|3T1) . It is different from m| 2 because for such a small eigenvalue the influence of the 
evolving background is relevant. The adiabaticity condition is still broken for m 2 ^^. The ef- 
fective mass squared is not always positive, but in the Bogolyubov coefficients method 
[321 1331 [34] that we use for calculating particle production the effective squared mass has to 
be positive only at times to and t\. For the previously discussed light mass matrix eigenvalue 
(|4T1) the adiabaticity parameter is smaller than for m| 2 (though it rapidly becomes of order 1 
during inflaton oscillations) and for all other eigenvalues of M' 2 (|3H the adiabaticity parameter 
is smaller than 10 -6 . This result means that the excitation £2 is the main channel for preheat- 
ing. This channel does not appear in global supersymmetry without non-renormalizable terms 
[221 ES], because in such a framework the eigenvalue related to the excitation around the phase of 
the flat direction is identically massless. The calculated numerically energy density of produced 
particles - e^ 2 grows very rapidly and begins to be comparable to the total energy density of 
all classical VEVs - ecLASS as can be seen in figure (fT6l) . The back-reaction of the produced 
particles on the classical VEVs evolution is not taken into account, so the result shown in figure 
(|16|) should be treated as an approximation only. However, from the approximate calculations it 
can be concluded that both non-perturbative particle production from the flat direction and pre- 
heating from the inflaton into particles corresponding to excitations £2 and (£7 + £12 + £15) /v3 
is efficient enough in order to melt the flat direction VEV, which unblocks all other channels 
of preheating from the inflaton. Moreover, to obtain this result one needs only the mechanism 
of particle production due to changing mass matrix eigenvalues related to the flat direction. 
The mechanism of particle production due to changing mass matrix eigenvectors [221 GS] is not 
necessary. 

Even though the phase dynamics is negligible in this scenario it is worth making one comment. 
In the diagonalizing base of excitations matrix A appearing in eq. (|28|) . has a non-trivial form 
only for excitations £2 and (£7 + £12 + £15) /v3 

. / cosa sina 
A=[ 

\—sina cosa 

With a = const = A = 1 and so M' 2 = M 2 , the mass matrix eigenvectors are constant in 
time and there is no preheating of the type proposed in [22] • However if the phase dynamic 
were present, the A matrix would mix in time two eigenstates with non-zero eigenvalues, which 
would lead to preheating from time-varying eigenvectors even from a single flat direction, in 
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3,1525xl0 7 3,155xl0 7 3,1575x 10 7 3,16x 10 7 
t 



Figure 16: Energy density of created particles e^ 2 in comparison to the total classical energy 
density of all VEVs ecLASS- Time on horizontal axes is expressed in Planck times. 

contradiction to [22] and [21]. This happens due to the inclusion of supergravity corrections or 
non-renormalizable terms. Both these effects make the scalar potential dependent on the flat 
direction, which gives masses to eigenstates £2 and (£7 + £12 + £15) /\/3- Without these effects, 
as in [22] and [25], eigenstates £2 and (£7 + £12 +£15) /V3 are massless and their mixing does 
not lead to particle production. 

4.2 X a = 1 and A x = 1 

In this scenario both udd and Higgs directions have non-zero VEVs during inflaton oscillations, 
breaking SU(3) x SU(2) x U(l) — > U(l). After eliminating Goldstone bosons in the unitary 
gauge we are left with the following parametrization of excitations 



\ (£3 + iU) e 




2 — 



5 (6 - *&) e 





Ui3 = \ (£l0 + e * CT ) 



(45) 



d j3 = I (£13+ ^14) e^, 
4 1 = 1 (£10 - «£li) e 4<T , 
d k 2 = \ (£13 - i$L4) e i<T , 
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Due to a more complicated set of excitations in this scenario the eigenvectors and eigenvalues 
of the mass matrix have been found numerically. The mass matrix is block-diagonal - one com- 
plicated block is related to excitations around non-zero VEVs and a separate, diagonal block is 
related to excitations £3, £4, £g, £9, £10, £11, £13 and £14 around VEVs equal to zero. 
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V 

Mass matrix eigenvalues corresponding to £3, £4, £s, £9, £10, £11, £13 and £14 are heavy as 
they correspond to Higgs particles related to breaking of non-diagonal generators. Their time 
evolution is strongly dominated by the udd flat direction or the Higgs direction VEVs - these 
eigenvalues are large and evolve slowly in time, which effectively blocks preheating of the inflaton 
into particles corresponding to these eigenvalues. The time evolution of some of the eigenvalues 
from the block corresponding to excitations around non-zero VEVs is also determined by large 
udd and H u Hd VEVs. The evolution of these eigenvalues does not allow either non-perturbative 
particle production from the flat direction or preheating from the inflaton (these eigenvalues 
correspond mainly to excitations related to Higgs particles of the diagonal generators breaking). 
This type of behavior was predicted by ref. [21]. Due to non-trivial phase dynamics the mass 
matrix eigenvectors also evolve in time allowing for non-perturbative particle production from 
flat directions as predicted in ref. [22J, which can lead to fast decay of flat direction VEVs. There 
is however another, more efficient channel of non-perturbative particle production due to the 
existence of light, non-adiabatically changing eigenvalues of the mass matrix. These light eigen- 
values appear due to the same mechanism as described in the previous case - they correspond to 
a combination of naturally light excitations around VEVs of complex fields a and x parameter- 
izing the (quasi) flat directions. In this scenario however, due to the non-zero VEV of the H U H^ 
direction, preheating from the inflaton is allowed from the beginning of inflaton oscillations into 
excitations around both directions. As an example figure (fT7l) shows the time evolution of such a 
light eigenvalue corresponding mainly to excitations around H u Hd direction. A comparison with 



m5~ 




Ml 



-5.X10 



Figure 17: Evolution of the light eigenvalue in comparison with the character of the evolution 
of the inflaton (red line) 



the time evolution of the inflaton VEV clearly shows that the time evolution of this eigenvalue 
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is dominated by the behavior of the infiaton. Moreover, because this eigenvalue is very small 
the impact of infiaton oscillations makes this eigenvalue periodically negative leading to a very 
effective tachyonic [35] preheating from the infiaton. 

5 Conclusions 

Achieving large flat direction VEVs through classical evolution during inflation is natural in a 
supergravity framework with non-minimal Kahler potential. Such large VEVs can block preheat- 
ing from the infiaton into certain channels. However supergravity effects and non-renormalizable 
terms, which create a potential for the flat direction, are a source of light, rapidly changing 
eigenvalues of the mass matrix. They allow the non-perturbative production of particles from 
the flat direction and preheating from the infiaton. Non-zero VEVs of Higgs fields can also lead to 
the existence of light, rapidly evolving eigenvalues, allowing for preheating from the infiaton into 
Higgs particles from the beginning of infiaton oscillations. Non-perturbative particle production 
due to the time evolution of the mass matrix eigenstates is not necessary to reduce flat direction 
VEV and unblock preheating. Thus non-perturbative particle production from the infiaton is 
likely to remain the source of preheating even in the initial presence of large flat direction VEVs. 
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